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Abstract. In the Gribov-Zwanziger scenario the confinement of gluons is attributed to an enhancement of 
the spectrum of the Faddeev- Popov operator near eigenvalue zero. This has been observed in functional and 
also in lattice calculations. The linear rise of the quark-anti-quark potential and thus quark confinement 
on the other hand seems to be connected to topological excitations. To investigate whether a connection 
exists between both aspects of confinement, the spectrum of the Faddeev-Popov operator in two topological 
background fields is determined analytically in SU(2) Yang-Mills theory. It is found that a single instanton, 
which is likely irrelevant to quark confinement, also sustains only few additional zero-modes. A center 
vortex, which is likely important to quark confinement, is found to contribute much more zero-modes, 
provided the vortex is of sufficient flux. Furthermore, the corresponding eigenstates in the vortex case 
satisfy one necessary condition for the confinement of quarks. 

PACS. 11.15.-q Gauge field theories - 12. 38. Aw QCD: Confinement and topological excitations 



1 Introduction 

The confinement of colored objects in QCD is still a chal- 
lenging problem, although much progress has been made 
in its understanding during the last few years QE]. This 
progress has been made along two directions, which seem 
to be quite different at first sight. One is a confinement 
scenario based on topological defects. The other is based 
on the properties of field configuration space, in the frame- 
work of the Gribov-Zwanziger scenario. It is yet unclear 
in which way both aspects are connected. The work pre- 
sented here is aimed at an investigation of a possible link, 
based on the following set of observations. 

The confinement of quarks is usually observed by the 
existence of an - up to string breaking - linearly rising 
quark-anti-quark potential. Already more than thirty years 
ago it was conjectured that topological excitations could 
be responsible for this behavior. Today, substantial evi- 
dence in favor of this scenario exists 0. This would be 
a very attractive realization, because such excitations can 
carry topological charge and hence be brought in contact 
with the spectrum of the Dirac operator via the Atiyah- 
Singer index theorem [3] and thus to chiral symmetry 
breaking by the Banks-Casher relation 0]. 

However, the nature of these excitations or correspond- 
ing topological defects is currently not finally resolved, nor 
is its gauge-dependence fully understood. Two of the most 
important candidates for such configurations are monopoles 
2,5 and center vortices 2,6,'Zj. Vortices seem to have 
slightly more attractive features in general [2], and hence 
will be investigated here. 



By construction, a non-interacting random ensemble of 
vortices yields an area law [Sj. Consequently, it has been 
shown by several groups independently that removal of 
center vortices removes quark confinement j^BH]. In addi- 
tion, this also restores chiral symmetry |10j . Thus, vortices 
will be used here as an example of a quark-confining topo- 
logical field configuration. 

On the other hand, instantons |11U12| . which play a 
role in chiral symmetry breaking and hadro-dynamics |12l 
IT5] , do likely not contribute directly to quark confinement 
and will be used here as an example of a quark non- 
confining topological excitation (but see e. g. ^3] on the 
topic of instantons and confinement). 

The confinement of gluons is to some extent a more 
subtle issue. As gluons carry adjoint color charges, string 
breaking is always present. Thus, it is much more com- 
plicated to judge whether gluons are confined or not. An 
empirical criterion for gluon confinement is the Oehme- 
Zimmermann super-convergence relation |15| . It states that 
a gluon is confined if its propagator vanishes at zero mo- 
mentum. Such a behavior has been observed in Landau 
gauge in calculations based on Dyson-Schwinger equations 
[Tlll6lll7| and renormalization group methods ^H] . Corre- 
sponding observations have also been made in Coulomb 
gauge using variational methods JU] . Lattice calculations 
have significant problems in reaching sufficiently large vol- 
umes (see e. g. [201 ) ; but there are indications on the 
largest, but highly asymmetric, lattices available [21], which 
support a vanishing gluon propagator. In three dimen- 
sions, where a similar behavior is found in Dyson-Schwinger 
calculations |17p22j . much larger lattices are possible. There, 
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more substantial evidence has been found in favor of such 
an infrared behavior of the gluon propagator |23j . 

This behavior is predicted by the Gribov-Zwanziger 
confinement scenario 17, 24 25 ,26 . The scenario starts 
from the Gribov problem 24,27]: Due to the gauge free- 
dom in QCD, not all possible field configurations are inde- 
pendent. It is therefore necessary to restrict the configu- 
ration space. Usual local conditions like Landau gauge are 
only sufficient in perturbation theory. Gribov proposed to 
restrict further to the first Gribov horizon, a compact re- 
gion around the origin, hence including perturbation the- 
ory. This is the domain in field space where the Faddeev- 
Popov operator 1 

M ab = -d„(d fi S ab +gf abc A c fi ) (1) 

is positive. At the boundary, a zero eigenvalue appears, 
which becomes negative outside the first Gribov region. 
Here g is the gauge coupling and f abc are the structure 
constants of the gauge group. This is not sufficient, as 
gauge copies exist inside this region |2B]- It is thus nec- 
essary to restrict even more to the fundamental modular 
region |29j. which has in part a common boundary with 
the first Gribov region. Up to now, no local characteriza- 
tion of the fundamental modular region exists. A non-local 
possibility |2S] for such a characterization is to choose the 
configuration on a gauge orbit which absolutely minimizes 

l#x{Al-A«^)(A%-Af*). (2) 

Herein AJJ is a n-instanton configuration to include fields 
in different instanton sectors. 

The Gribov-Zwanziger scenario utilizes the entropy ar- 
gument that the volume of an infinite dimensional volume 
is concentrated at its border. Thus, the configurations at 
the common boundary should dominate the partition sum 
|25| . Hence, configurations with a vanishing determinant 
of the Faddeev-Popov operator dominate the infrared, 
and the spectrum of this operator should have an enhance- 
ment of zero-modes or near-zero modes compared to the 
vacuum case. 

In Landau gauge, the expectation value of the Faddeev- 
Popov operator is the inverse ghost propagator. Thus, the 
scenario predicts an infrared enhanced ghost propagator. 
It can be shown quite generally that this prediction is 
correct |3()| : it has also been confirmed in many calcu- 
lations [Toll 1 711 18ll2l)] . including lattice ones |3*T1I^] . The 
connection between such an enhancement of the spectrum 
and an infrared divergent ghost propagator has also been 
shown explicitly in lattice calculations |33|. Furthermore, 
the scenario predicts that the infrared behavior of Yang- 
Mills theory is dominated by the gauge-fixing part of the 
action |25|. which has been confirmed in Dyson-Schwinger 
calculations ^H]- In addition, the properties of the ghost- 
gluon vertex have been found to agree very well with this 
scenario |25II34| . Similar relations hold in Coulomb gauge, 

1 As the topological field configurations are defined in Eu- 
clidean space, it will be used exclusively here. 



which have also been confirmed in functional and lattice 
calculations |19U35| . 

Corresponding predictions are also given by a scenario 
proposed by Kugo and Ojima |3ri) . which is based on the 
BRST symmetry. Especially, at least in Landau gauge, an 
enhancement of the ghost propagator is predicted |37J , and 
therefore an enhancement of additional zero-modes in the 
spectrum of the Faddeev-Popov operator. It is, however, 
not yet clear, if a connection between both scenarios exist. 

These observations let it appear very promising to in- 
vestigate analytically which type of connection between 
quark-confining topological field configurations and an en- 
hancement of zero or near-zero modes of the Faddeev- 
Popov operator exists. This is also motivated by indepen- 
dent results on this topic in lattice gauge theory [123 . As a 
cross-check, it is interesting to study a configuration like 
an instanton, which is probably not directly involved in 
quark confinement. 

This will be the content of this work: In section the 
spectrum of the Faddeev-Popov operator in the vacuum 
will be discussed, in order to fix notations, for compari- 
son, and for the sake of completeness. The spectrum in 
a one-instanton background will be determined in section 
13 and for in a one-center- vortex background in section 0] 
In section |5] a necessary criterion for confinement will be 
tested for the solution in the vortex case. The results and 
some possible implications will be discussed in section 

2 Vacuum 

In the vacuum the Faddeev-Popov operator 0) is just the 
negative of the 4-dimcnsional Laplacian times a unit ma- 
trix in color space, 

/ n2 n2 n2 n2 \ 

= S ab (Xrdydr + - 1 g d„ Sin 2 V dr, 

V r J r 2 sin r\ 

+ ^^(-^dl + ^-d e smede)) (3) 
r 2 sur T] V sin 6 v sm 6 J J 

= -S ab f-d r rd r + \dl + -d p pd p + i# ) . (4) 

The second expression is given in hyper-spherical coordi- 
nates 

r fl = (x,y,z,t) T (5) 
= (r cos 4> sin 9 sin r], r sin <fi sin 9 sin 77, 
r cos 9 sin 77, r cos r]) J ', 

where <fi ranges in [0, 27r) and 6 and rj range in [0, tt). The 
third expression is in bi-polar coordinates, essentially a set 
of two polar coordinate systems, 

77, = (r cos9,r sin9, pcosr], psmrj) T . (6) 

In the vacuum the eigenmode equation is separable and 
thus directly solvable. 
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In Cartesian coordinates, the eigenmodes, solutions of algebra elements by j3H] 
the eigenequation 

A u = - 



2 m 



M a0 4>° = oj z 4>' 

of the vacuum Faddeev-Popov operator are plane waves 
C a exp(iA; /J r Al ), and the spectrum is the continuum of pos- 
itive real numbers ui 2 with ui 2 — k^k^. Note that in the 
vacuum the operator is positive definite. As these are free 
states, they are non-normalizable and correspond to scat- 
tering states in quantum mechanics. In the vacuum the 
colors always decouple, and there is always a set of N 2 — 1 
independent solutions. The only zero-modes are trivial, 
constant functions, of which also N 2 — 1 linearly indepen- 
dent ones exist. 

In hyper-spherical coordinates, the eigenmodes can be 
decomposed into partial waves with three independent in- 
teger angular quantum numbers n, I, and m |38| . The 
eigenmodes are then given by 



<t> a (r» 



C 



(cos 2 77 — 1) 



1/4 



Y' 



where J are the Bessel functions and P are the associated 
Legendre functions of the first kind. Y^ are the ordinary 
spherical harmonics. A non- vanishing zero- mode is only 
possible for n = I = m = 0. These are again constant 
solutions. 

However, further zero-modes, e. g. of the form C a /r 2 , 
also exist. These will not be admitted as they are in gen- 
eral more singular than the gauge field configuration. This 
will also be required when treating non-zero gauge field 
configurations below. This will imply that in the cases 
treated here any admissible solution may not diverge at 
spatial infinity. 

In bi-polar coordinates, the problem separates into two 
two-dimensional sets, one for each polar coordinate sys- 
tem. As both angular coordinates are 27r-periodic, the so- 
lution is 



C a J\ 



(^uj \ZT~~s 2 p^j J| m |(wsr) exp(i(m9 + nrj)).(7) 



J n are again the Bessel functions. The angular quantum 
numbers m and n are positive and negative integers in- 
cluding zero and can be chosen independently for each 
color. The continuous variable s, which ranges in [0,1], 
"shifts the eigenvalue" between both coordinate sets. Again 
there are only trivial zero-modes, all constants: For u> = 0, 
only n = m = leads to non- vanishing values of the Bessel 
functions. 



3 Instanton 

3.1 Field configuration 

The simplest case of a topological field configuration in 
Yang-Mills theory is an instanton 2 . These are given as 

2 The remaining part of the article is restricted to SU(2) 
Yang-Mills theory. 



A 2 



' \1V ' V 



T M = (lT, 1) 
V = (-*T, 1), 

where A characterizes the size of the instanton and t 1 are 
the Pauli matrices. The corresponding gauge fields are 
given by the regular functions 



AO. 



1 



g r 2 + A' 



~ r vCvu 



(8) 



The constant real matrices £° (the 't Hooft tensors) form 
the algebra 



{C a ,C b } = -5 ab - 



(9) 
(10) 



The instanton fields are transverse, d^A^ = 0, and 
thus admissible in Landau gauge. Furthermore, they are 
part of the fundamental modular region as they trivially 
minimize (j2Jl. Hence, the Faddeev-Popov operator should 
have only positive or zero eigenvalues, which will be con- 
firmed below. 



3.2 Analytical treatment 

Due to the structure of (JSJ and the transversality of the 
field, the eigenvalue equation of the Faddeev-Popov oper- 
ator in an instanton field can be written as 



8 



2 j,a 



f 



abc 



A 2 



2ia 



(ii) 



Neither do the different colors decouple nor is the prob- 
lem fully separable anymore. As it is possible to express 
distances r in dimensionless variables r/A, only the dimen- 
sionless quantity ui\ is left as an independent parameter. 

First of all, due to the transversality of the instanton 
field, the three trivial constant zero-modes still exist. 

To solve the set of equations (fTTJ) . the first step is to 
note that the expressions r^^du are angular momentum 
operators. The eigenvalue equation can then be rewritten 



= (d 2 + uj 2 )^ 1 
2 



(12) 



. , . 2 A{L M + L 12 )4> 2 + {-L 2i + L 13 )(j) 3 ) 
z(r z + A z ) 



= (d 2 +cj 2 )^ 2 
2 



(13) 



+ — —(-(L34 + L 12 )4> x + (L u + L 23 )(/) 3 ) 

(14) 



i(r 2 + A 2 ) 

= (d 2 + cj 2 )<P 3 
2 



+ -, 



i(r 2 + A 2 ) 
L a p = i(x a d/3 - xpd a ) 



((L24 - - (L u + L 23 )cj) 2 ) 
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The angular momentum operators satisfy the commuta- 
tion relations 

[Lap, Ljg] — i(8p*,L a 5 + SasLpy + 5psL ia + S aj Lsp). 

This implies that any solution which only depends on r 
does not feel the presence of the instanton, and is the 
same solution as in the vacuum. Thus, there exists a con- 
tinuous positive spectrum, as expected. Furthermore, the 
Laplacian can be written as |38| 



d 2 = ^d r r 3 d r - i 



l<a</3<4 



(15) 



As only three particular linear combinations of the six an- 
gular momentum operators appear, it is possible to define 
the set of operators 

L 1 = — (L14 + L23) 

L 2 = — (£13 — L24) 

L 3 = - (£34 + £12), 
which form a three-dimensional spin-algebra 

[L a ,L b ] =if abc L c . 
In addition, explicit calculation shows that 
L 2 = (L 1 ) 2 + (L 2 ) 2 + (L 3 ) 2 

l<a</3<4 

+ 2 {L14L23 — L13L24 + L12L34). 

The second term vanishes. Thus, the angular part of the 
Laplacian operator can be replaced by AL 2 . 

d 2 = ^d r r 3 d r - 

Hence, the problem can be separated in a radial part and 
an angular part. The latter can be characterized by two 
independent quantum numbers I and m. Selecting L 3 as 
"quantization" axis 3 , the eigenfunctions depend on /, de- 
noting an "orbital angular momentum" , which can be in- 
teger or half-integer. The "magnetic" eigenstates are then 
counted by m. In principle, each color could have its own 
angular momentum /, but in the equations (|12I14|) no op- 
erator appears which can change /. So to compensate the 
angular dependencies, each color has to have the same l- 
value for a solution, and L 2 can be replaced by /(/ + 1) as 
its eigenvalue. 

3 This is the axis characterizing rotations in the 1-2 and 3- 
4 planes around the same angle. In fact bi-polar coordinates 
would be the natural way to solve this problem explicitly in 
terms of r 2 + p 2 , r/p and 9 + rj. However, this will be not 
necessary for the purpose at hand. 



Thereby each of the equations splits in two parts, 



(L 3 6 2 +L 2 6 3 ) 



+ 



i{r 2 + X 2 ) 



, 1 „ 4/(7 + 1) ., , ., 



+- 



i(r 2 + A 2 ) 



-L 



3j,1 1 rl^3N 



+ 



i(r 2 + X 2 ) 



<L 2 tf - LU 2 ) 



Using a matrix representation for the angular momentum 
operators, this system can be written as a matrix equation 



(lD r + L I )cp = 0. 



(16) 



The dimensionality of the vector is 3(2/ + 1), repre- 
senting the three colors and the (21 + 1) independent m- 
quantum numbers. D r collects the radial derivatives and 
terms after multiplication with r 2 + A 2 . 1 denotes the unit 
matrix in this space and Lj is the hermitian matrix 



(17) 



Expanding the solution cp in the eigenbasis {li} of L T as 
^2<j)ii(r)li completely decouples the angular and radial 
parts, leading to the radial equation 








L 2 


L 3 





L 1 


L 2 


-L 1 






(D r + Ci)4> u (r) = 0. 



(18) 



Here, Cj are the eigenvalue of the corresponding eigenvec- 
tor. The eigenvalue problem is a typical spin problem. In 
general, there are 3(2Z + 1) eigenvalues, which can be either 
positive or negative or eventually zero. Note that in the 
case / = the solution are functions dependent only on r, 
as already discussed above. Therefore, only Z > 1/2 will 
be regarded here. The eigenvectors will not be needed to 
determine the spectrum, but could be obtained in a direct 
way. However, to solve the radial equation i|18|) . the eigen- 
values have to be determined. A close inspection of the 
eigenvalue problem yields the characteristic polynomial of 
(fTTjl as 



(c - 4) 2Z+1 (c + 4Z) 2i+3 (c - 4(1 + l)) 5 



= 0. 



Hence there are only three different eigenvalues c^: 4, —4/, 
and 4(1 + 1) with multiplicities 21 + 1,2/ + 3, and 21 - 1, 
respectively. The radial solution does not depend on the 
multiplicity index and <pu can be replaced by <pi c . With 
this, the radial equation can finally be written down as 



= \d r r 3 d r (j> lc +[lu 2 - 



41(1 + 1) 



A 2 



<M19) 
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The only difference between the different colors are then 
their components along the quantization axis and thus 
their composition in terms of the different magnetic eigen- 
states: The different colors are rotated with respect to each 
other. 

A first step to solve the ordinary differential equation 
for the radial part 1)19(1 is to investigate the asymptotic 
properties in order to study the existence of admissible 
solutions. At small r, the only relevant part is the angular 
part, leading to the simpler equation 



= d r <f>i c + -d r <j>i c 



4/(1 + 1) 



and its convergence 4 cannot be checked easily for r > A. 
It is therefore necessary to perform a numerical analysis, 
which will be done in section l3~3l 



3.3 Results 

Beforehand, it is worthwhile to treat the case uj 2 = 
explicitly. The equation 



A 2 



This equation is independent of c and uj 2 . It is solved by 

with C\ and Ci integration constants. Therefore, there ex- 
ists only one solution which is regular at the origin, and 
there is no more than one non-singular solution for the dif- 
ferential equation l|19|) . The homogeneity of the differential 
equation then implies that the only remaining integration 
constant is an over-all normalization factor. 

At large distances only the part containing the eigen- 
value in 1)19(1 is relevant as long as to 2 is non-vanishing. 
Thus, the more simple differential equation 



d r (f>ic + -d r 4>i c + uj 4>i c 
r 







has to be treated. This is just the equation for 1 = 0, which 
is solved by Bessel functions of the first and the second 
kind. In this case the second solution is also admissible, 
as the behavior at the origin is not relevant, yielding 



D Uur) +D Y n {u:r) 



(20) 



D\ and Z?2 are new integration constants. In principle, 
1(20(1 allows to construct finite solutions for any positive or 
negative uj 2 by an appropriate choice of D\ and D^. Thus, 
besides showing that the long-range regular solutions will 
be of type cos(r/A + <5)/r 3 / 2 , 1(20(1 does not provide fur- 
ther information. It is therefore necessary to obtain an 
explicit solution for the differential equation ((19(1 . Making 
the ansatz <\>i c = r 2l wi c (r), and expanding then wi c (r) in 
powers of r/X, this is straightforward. The general solution 
for any w 2 , I, and c is given by 



kc = r 21 o„ 

n=0 



r 
A 2 



2 \ n 



(21) 



a_x = 
a = D 

_ u 2 A 2 (q n _i - a n _ 2 ) 
fl " ~~ 4n(n + 1) + 8ln 

t ((4n + 8Q(rc- 1) + c)a n _i 
4n(n + 1) + 8ln 

Here D is the free overall normalization factor. Unfortu- 
nately, this series is similar to the hyper-geometric series, 



is directly solvable, and yields hyper-geometric functions 

v2\ 1+* 



(r)=C 1 (jz) 2 F 1 (-i-/-i v / l-c + 4Z + 4Z 2 , 
-^-l + ^Vl-c + 4l + 4l 2 ,-2l,-^ 

„2N I 



c + Al + Al 2 , 



1 1 / r 2 

- + / + -V1-C + 4/ + 4/ 2 , 2 + 21, - — 

2 2 A z 



As 2 fi(a, /3, 7, 0) = 1, the first of these solutions is singu- 
lar at the origin and can be dismissed. The second solution 
is regular. Thus, it remains to investigate its long distance 
behavior. Using an identity for the hyper-geometric func- 
tions |40|. the expression can be rewritten as 



r 

A 2 



-(i-H-iVl-c+4!+4Z 2 ) 



X 2 

3 

2 " 



F^+l-^l-c + Al + Al 2 , 
f-J- -Vl -c + 4Z + 4P,2 + 2/,^2 

2 £1-1 



For r — > oo, the hyper-geometric function is then a con- 
stant for all eigenvalues. The pre- factor is finite at spatial 
infinity only for I = 1 and c — 8 with multiplicity one and 
I = 1/2 and c = 4 with multiplicity two. In these cases 
the zero-modes can be given in closed form and read for 
1 = 1/2 



4> hA {r) = 2D- 
and for I = 1 

0i 8 (r)=3U- 



-h+ 1 + In 1 + ^r 



-2(l + ^)ln 1+?, 



A 1 



4 To be precise: As this is a solution to an initial value prob- 
lem for a smooth differential equation, the series is guaranteed 
to exist and converge for any finite r. This, however, does not 
prevent a divergence for r — > oo. 
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No further zero-modes exist for larger I or for any other 
eigenvalue c. Hence, the instanton sustains six zero-modes 
in total, the three trivial ones at I — 0, two for 1 = 1/2 
and one for I = 1. 

These results can be supported by numerical calcula- 
tions. A selection of numerical results 5 is shown in figure^ 
for 1 = 1/2 and in figure [21 for I = 1 for the possible values 
of c and various values of uj 2 X 2 . Calculations at larger / 
or other values of ui 2 \ 2 confirm these findings. No further 
zero-modes exist. 



3.4 Singular gauge 

It is an interesting question, to which extent these findings 
are gauge dependent. In general, different gauges lead to 
much more complicated differential equations, especially if 
the field configuration is no longer transverse afterwards. 
In the case of the instanton, it is possible to check at least 
one other gauge. Performing the gauge transformation ^Hj 

G(x) = 

r 

the instanton field configuration is given by |39j 
- 2A 2 



r 2 (r 2 + A 2 ) 

This gauge transformation has the property to move the 
non-trivial behavior from infinity to a finite region, at 
the cost of a small-distance divergence. In a sense, this 
gauge transformation is one of the most severe changes 
which can be done to the global, and thus long-range and 
confinement-relevant, properties of a gauge field configu- 
ration. Therefore, it is an important test of the results 
found above. 

The problem can be treated along the same lines as be- 
fore, but as the field configuration now diverges as 1/r 2 at 
the origin, corresponding eigenfunctions have to be admit- 
ted. It is then still possible to separate the angular struc- 
ture, and the equations turn out to be the same as (jlfjjl . 
except for a different quantization axis (which nonethe- 
less can be called L 3 ). In addition, an additional factor of 
r 2 /A 2 appears in the radial part. Thus, the same angu- 
lar solutions are found, and the radial equation takes the 
form 







1 



Die 



41(1+1) 



cA 2 



r 2 (r 2 + A 2 ) 



•He 



5 These were obtained by summing up the series explicitly 
for r / A < 1 . At some point r s , the series was used to find initial 
values for the differential equation, which was then solved by 
a Cash-Karp Runge-Kutta algorithm |41| . This was necessary 
due to explicit divergencies of the angular term in the differ- 
ential equation at 0. Otherwise the numerical solution could 
have been obtained by the integration method alone. 



Radial eigenfunction at 1=1/2 and c=-4l j 





--crA 2 =-1 






-co 2 ^ 2 =0 
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Radial eigenfunction at 1=1/2 and c=4 j 
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Fig. 1. The radial eigenfunctions <f>i c for 1 = 1/2 and the two 
different c values. For better visualization, positive w 2 -solutions 
have been normalized so that their maximum is 3, while modes 
with u> 2 \ 2 < have been normalized so that (f> lc /r 2l \,— o = 1. 



with the same values of c as before. The solutions of this 
equation differ qualitatively from their previous counter- 
parts only by an admitted divergence at zero, and are 
given in explicit form for u> 2 ^ by 



= -l+V(2'+l) 2 -c 



71=0 



a_i = 
a Q = D 



(((c-\W)-4(n-l) 2 

-2(n-l)V(2; + l) 2 -c)a„_! - A 2 w 2 a„_ 2 ) 



2n (2n + ^(21 + l) 2 - c) 
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| Radial eigenfunction at 1=1 and c=4(l+1) | 
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| Radial eigenfunction at 1=1 and c=-4l | 
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Radial eigenfunction at 1=1 and c=4 | 




Fig. 2. The radial eigenfunctions (f>i c for I = 1 and the three 
different c values. For better visualization, positive oj 2 -solutions 
have been normalized so that their maximum is 3, while modes 
with u> 2 \ 2 < have been normalized so that (f>i c /r 2l \ r =o = !• 



There are three non-trivial zero-modes, at exactly the 
same values of I and c, reading 



^l+(l + ^)ln(^) 



1 + 2^ + 2 



3 4 2 
k 8 = D 



These solutions decay faster than in the regular case, and 
are even square- integrable and therefore normalizable and 
localized. 

Note that there are now three additional zero-modes 
at Z = 0, <p(r)o o = C a /r 2 , which previously have been 
rejected in a non-singular field-configuration. These arc 
also trivial zero-modes, in the sense that they already exist 
in the vacuum, and the number of non-trivial zero-modes 
is the same. 

Hence, the non-trivial properties of the spectrum is 
not affected by this gauge transformation. This is not en- 
tirely expected, as the Faddeev-Popov operator is gauge- 
dependent, and the chosen gauge transformation severely 
affects the long distance properties of the field configura- 
tion. 



4 Center Vortex 



4.1 Field configuration 



In this section the spectrum of the Faddeev-Popov opera- 
tor will be studied in the background of a thick, oriented 
center vortex. The vortex field used here provides a non- 
trivial Wilson loop in the fundamental representation of 
the gauge group. The field strength of such a vortex in 
SU(2) is given by @21 



A% = 6 



3a 1 P(p) 
9 P 



All other components vanish 6 . This field is transverse. The 
function jj,(p) varies from zero at p — to zero at p = oo 
for a flux zero vortex. For higher fluxes, it varies from zero 
at p = to the flux 2n + 1 at p = oo. Herein, n is a positive 
integer or zero. There is no further specification of this 
"profile" of the vortex. For the results, several different 
"profiles" will be compared. The only assumption made 
will be that they are smooth. 

This type of center vortices are energy minimizing con- 
figurations after including first order quantum effects |42| . 
At the classical level, only instantons are such solutions. 



In principle, it would be possible that there is a non- 
vanishing component A p . Such a component can be removed 
by choosing a suitable gauge |42) . 
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4.2 Analytical treatment 

The natural coordinate system for this problem are bi- 
polar coordinates JBJ. The three eigenvalue equations for 
the three colors are then given by 



aV - ^d v <t> 2 = W V 



ii 
p 2 ' 
p_ 



-9V + iia/ = w ^ 2 
p 2 

-d 2 (/) 3 = uj 2 4> 3 . 



The equation for the 3-component decouples and thus is 
solved by the free solution (7J). This implies that </> 3 be- 
comes constant for any zero-mode. It thus remains to solve 
the coupled equations for the eigenfunctions of colors one 
and two. For these the behavior in r — 9, in which they 
are trivial, and p — rj decouple. The equations are then 
reduced to a two-dimensional problem 

(jd p pd„ + j 2 dl + u 2 (l - S 2 )) 1 + ^d^ 2 = 

\^d p pd p + ^d 2 +uj 2 (1- s 2 )) 2 - ^d^ 1 = 0. 

As the angle i] is 27r-periodic, it is possible to expand the 
eigenfunctions in a Fourier-series as 4> a — J2 m c m exp(imrj). 
This yields 

° = E(( 1 ^p- ! ^+- 2 ( 1 - s2 )) c ™ 



(22) 



imp, 2i, 
— — c m J exp(«m?7) = 



= E((^-^+- 2 d- 52 )) c 

2 " C m) exp(lTO77). 



imp 



(23) 



Only the equations for the same m are coupled, as r\ is 
arbitrary, and it suffices to solve the equations for each m 
separately. 

At to = the system decouples, as it is reduc ed to th e 
free case. Consequently, this yields eg = C° Jo (ojyjl — s 2 p) 
as a solution. This implies that s is again constrained to 
< s < 1, which is confirmed by the large p solution. 
Thus, to can be taken non-zero in the following. 

Since p goes to zero at the origin, the angular term 
to 2 //? 2 dominates, just as in the case of the instanton. 
The solution is then similarly given for small 7 p by 

c a m ~ A a p\™\ +B a p~^, (24) 

with integration constants A a and B a . As for all profiles 
the field configurations vanish at the origin faster than 

7 Small is here determined by the typical scale of the profile 
function p. Whenever there is a notion of large and small, 
this will be with respect to this intrinsic scale, which is only 
determined if the function p is given explicitly. 



1/p, the second solution is dismissed. This requires two 
of the four free (complex) integration constants. As the 
differential equations are homogenous in the functions, the 
remaining integration constants are again scale factors. 

At very large distances for non-zero w, the eigenvalue 
term dominates. The equations decouple and in this case 
the eigenfunctions will behave, as in the free case, like 
Bessel-functions. Thus, for finite uj the influence of the 
vortex is only relevant at intermediate distances. The sit- 
uation is different in the case of a zero-mode, and will be 
treated below. 

To make progress with the solution, it is useful to 
rewrite equations Ij22l23|l with explicit real and imaginary 
parts. Then the equations for the four independent real 
functions b'L + ie'L = dL become 



-dppdp - 
P 


9 

m 






1^- 


mp 2 

2 m 

p 


= (25) 


-dppdp - 
P 


o 

m 

T 2 " 


f^ 2 (i- 




1 eln + 


mp 2 

p 2 m 


= (26) 


-dppdp - 
P 


2 

TO 

T 2 " 






\ b 2 m + 


mp 

2 

P 


= (27) 


-dppdp - 
P 


2 

m 

I 2 " 






\e 2 m - 


mp ! 

2 m 

p 


= 0. (28) 



Highly advantageously, the equations decouple into two 
sets of equations, which only differ by the sign of to. There- 
fore b x m ~ b 2 _ m and e\ n ~ e 2 _ m , and it suffices to solve the 
first and fourth equation of 1)25128(1 . As now the last term 
has the same sign, it is again possible to decouple the 
equations by looking for solutions of type 6+ = 6^ + e 2 n 



and b m = b}„ — < 

by 



The original functions can be recovered 



Note that the scale constant for the two solutions b^ can 
be selected independently: It is possible to have b] n = 
ie 2 ^. The decoupled equations are then 



-dppdp -™+a J 2 (l-s 2 ))b+- ^bt = 



-dppdp - + c^ 2 (l - s 2 ) ) b~ + ^b~ = 0. 



These equations are the same up to the sign of to, yield- 
ing 6+ = bZ m . Note especially that for b x m and e 2 ^ to be 
convergent simultaneously, at least 6+ or 6~ has to be 
convergent. A cancellation of divergencies if both are di- 
vergent is not possible. Hence it finally suffices to solve 
the single equation 



dppdp 



in 



+ u 2 (l~s 2 ))b 



mp 

„2 " m 



b+ = 0. (29) 



The small- p solution for &+ is the same as before, (|24|) . 
At large p the same applies as previously, but it is now 
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bt = 0. 



more interesting to treat the zero-modes. In this case the 
equation at large p is given by (/i(oo) = 2n + 1 for a vortex 
offlux2n+l) 

'1 Q _ m 2 +m{2n+iy 
-d p pd p 

Since to is integer, there are three possibilities. 

First, if m 2 + (2n+l)m > 0, which is the case for m > 

or to < — (2n+l), the solution is given by p ± \ /m2 + m ( 2n + 1 ) , 
It then remains to determine whether the small-p con- 
vergent solution connects to the large-p convergent solu- 
tion or not. It will turn out that the convergent small- 
p-solution connects to the divergent large-p solution and 
vice versa. Thus this type of solution cannot lead to an 
admissible zero-mode. 



Second, if (2n + l)m 



which is possible for neg- 



ative, odd to, the second term vanishes, and the solution 
is of type c + ln(r). It is again important to determine 
whether the small-p convergent solution connects to the 
constant or divergent solution. Again it will turn out that 
a small- p convergent solution connects to a large-/? diver- 
gent one. 

The last possibility is to 2 + (2n + 1)to < 0, which is 
the case for to negative and |to| < (2n +1). This is only 
possible for n > 0, and only flux three and higher vortices 
can produce such solutions. This yields a large-p behavior 
of 



bt 



C-y cos ( y/\m 2 + m(2n + X)\bx(p) 



+C 2 sin (V |to 2 + w(2n + 1)| ln(pfj 



Although this is not vanishing, it is finite. In these cases 
there exist finite zero-modes, as will be confirmed below. 
So, 2n permitted absolute values of m exist for each n. Due 
to the different angular structure, positive and negative re- 
solutions are independent. Hence, all in all 4n zero-modes 
for a vortex of flux 2n+l exist. E. g. for n = 1 there are for 
the +- function one for —1 and one for —2. For to = —3, 
to 2 + (2n + 1)to = 0. In addition, there are two for the 
— function for to = 1 and to = 2. Note that for n = 
the vortex has a non-vanishing flux of one but does not 
support any zero- modes. The angular quantum number of 
the r — r\ part of the solution has to be zero and thus does 
not increase the multiplicity. 

Furthermore, neither the small nor the long distance 
behavior depends on the the specific form of p(p): The 
existence and number of zero-modes is independent of the 
vortex profile, only the total flux is relevant. 

It is finally necessary to obtain the intermediate range. 
To perform this task, equation (|2~9T has to be solved. As 
in the case of the instantons, this is possible using a series 
expansion, provided a series expansion of p = pip 1 
exists (which will necessarily have po — 0). This condition 
is very reasonable for a smooth function going form to 
2n + l without any singularities. The solution of equation 
Ij29(l is then given by 



bt 



p |to| $>+„p" 

n=0 



(30) 



m—1 



b+ = 



-uj 2 {1 - s 2 )b+ n _ 2 + mJ2^ b+ m _ 1 _ i p 1+ 
n 2 + 2n|m| 



Unfortunately, as in the instanton case, this is again a 
series of hyper-geometric type, and it is not practical to 
use it for explicit calculations for large p. This will again 
be remedied by use of a numerical treatment. 



4.3 Results 

The calculations above are confirmed by a numerical cal- 
culation, using the same techniques as in the instanton 
case. A flux zero vortex of profile p — pX/(X 2 + p 2 ) and 
a flux one vortex of profile p = p/(X + p) do not sus- 
tain zero-modes. The results for a flux three vortex and 
a representative selection of to- and lo 2 (1 — s 2 )- values are 
shown in figure |3| It is clearly visible how at negative m 
the different convergent solutions appear, confirming the 
analytical calculations. This is even more evident in fig- 
ure 21 where the zero-modes for flux three and flux five 
vortices with the profile (2n + l)p/ (A + p) are shown as a 
function of r and to. 

Further calculations have been performed with other 
vortex profiles 8 , taken from [42] . e.g. p = exp(— p/A— X/p), 
p = (2n+l) exp(-A 3 /p 3 ) and p = p 6 /{X 6 +p 6 ). These cal- 
culations do not yield qualitatively different results than 
the above findings. 



5 A confinement criterion 

In it has been argued that the criterion 



lim — '-^ — '- > 



(31) 



is in Coulomb gauge a necessary criterion for confinement, 
although not sufficient. A fulfillment leads to a divergent 
energy of an unscreened color charge. This is of course not 
the ground state of QCD, but this must occur in a confin- 
ing theory. The quantity 7 is the density of eigenmodes at 
eigenvalue w 2 , and the function F is the expectation value 
of the negative Laplace operator in the eigenstates of the 
Faddeev-Popov operator to eigenvalue u 2 . The derivation 
of this condition only applies in Coulomb gauge. Thus, 
this criterion cannot be applied to the instanton configu- 
ration, but only to the center-vortex configuration, when 
the vortex is taken to be not aligned in the time direction. 

In the following the criterion will be studied for the 
vacuum, the instanton (which will be included for com- 
pleteness), and the vortex. 

8 Some of these profiles do not have a simple series expan- 
sion, so (I3UB is only of a limited value. Nonetheless, these are 
amendable to numerical calculations, and the results do not 
differ qualitatively from profiles where a simple expansion ex- 
ists. 
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| Radial eigenfunction at m=2 | 



Zero-modes for different m for a flux 3 vortex 
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Fig. 3. The radial eigenfunctions for different m and u (1 — 
s 2 )-values in a flux 3 vortex. For better visualization, positive 
w 2 -solutions have been normalized so that their maximum is 
3, while modes with uj 2 < have been normalized so that 
&+/p M |„=o = l. 
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Zero-modes for different m for a flux 5 vortex 




Fig. 4. The radial zero-modes fo+ for different m in a flux 3 vor- 
tex in the top panel and in a flux 5 vortex in the bottom panel. 
The solutions have been normalized so that b+/p' m '| p=0 = 1. 



5.1 Vacuum 

The eigenfunctions in the vacuum are plane waves. There 
are no non-trivial zero- modes, 7(0) = (in fact, it vanishes 
as ~ uj |35p. Hence, it remains to determine the function 
F. In 35 it is given by 

F(uj 2 ) = -W rf^C n(*)d 2 ^ „(*), 

n 

where „ is the n-th eigenstate to the eigenvalue uj 2 . In 
an infinite volume, this has to be normalized differently 

as 



F(uJ 2 ) 



The corresponding calculations for the eigenfunctions 
C° exp(ifc^x M ), with k 2 = uj 2 , can then be done explicitly 
to yield 



F{UJ) - EM 2 jd^jd^ ~ k 



(32) 
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Thus, F(uj 2 ) = uj 2 , as already stated in [jSJ Consequently, 
the condition ((31(1 is not fulfilled, as ^(uj 2 ) • uj 2 /ui 2 — as 
uj 2 -► 0, since 7 (0) = 0. 



5.2 Instanton 

In cases with a non-trivial gauge field, the situation is more 
complicated. In both cases treated here the eigenf unctions 
are eigenstates of the angular momentum operator, and it 
would be possible (and in fact has been done) to calcu- 
late the condition explicitly, at least for uj 2 = 0. This is 
not very enlightening, and a different way is pursued here. 
Note first that 7(0) 7^ 0, as there exist for the instanton 
as well as for the vortex non-trivial zero-modes. Further- 
more, as the states are eigenstates of the Faddeev-Popov 
operator, it follows that (for transverse configurations) 

- c* „ = u 2 <t>i% ji* „ + r bc A>>rj „■ 



F is then given by 



F(uj 2 



■ UJ 



y f <Px f abc A b n 6 a l d u 6 c 2 



uj' 1 n ^ ui z n 



and the vacuum result follows directly without any calcu- 
lation. It is furthermore clear that if the gauge field and 
thus the eigenmodes are regular everywhere, the second 
term has to be behave as —uj 2 for ui 2 — > in order to 
prevent the fulfillment of ((31J) in case of a non-zero level- 
density at uj 2 = 0. 

In the instanton case, the eigenmodes are eigenfunc- 
tions of f abc (x 2 + \ 2 )A b d a . Using then the orthogonality 
of the angular part of the eigenfunctions to perform the 
angular integrations yields 



F(uj 2 



UJ 



j: n fx 3 dxr 



be 4c n 



X 2 +U! i 



uj* n~ui* 



En / X*dx4> a 



uj* n T ui* 



With the known exact formulas for the radial part, the 
integrations can be performed exactly, but yield a lengthy 
expression involving several poly-logarithms. This expres- 
sion, however, vanishes in the infinite volume limit, in ac- 
cordance with a numerical integration. In this case F{ui 2 ) = 
uj 2 , as in the vacuum. Although now 1)31(1 is fulfilled, noth- 
ing is implied by this, as the instanton is not a Coulomb 
gauge configuration. Still, this result is not quite the one 
which would be naively expected for an instanton, but it 
is not so suprising when considering the presence of zero- 
modes. 



5.3 Vortices 

Different Fourier modes are not coupled by the gauge field. 
Therefore, it follows that 



/ 



abc^Pip) ,a* 



= .f 



a3cP(p) , a* 



u ui 2 mn u V'rui 2 mn 



im[ip , 2 * ,1 _ ,1* ,2 \ 

2 V^ui 2 mnrui 2 mn Yu) 2 nm^ui 2 mn) 



2mn{p) , &1 



e 2 , -b 2 ? e 1 , 

2 V u) van uj z mn to 1 * mn u)^ mn} 



In the last line the notation of the previous calculations 
has been used. Such contributions exist for all allowed m 
values. As the eigenfunctions are independent, but have a 
common long range behavior, the fraction in F only van- 
ishes if the denominator scales with a larger power of the 
volume than the numerator. Otherwise any desired result 
could be obtained by an appropriate choice of the inde- 
pendent normalization constants. In the denominator this 
is of no importance, as the integrated functions are posi- 
tive semidefinite. Thus, no singularity can be constructed. 
In addition, from the trivial two dimensions, each term is 
multiplied by a factor stemming from a free- wave solution, 
i. e. a term proportional to the square root of the volume, 
which cancels out. 

The zero-modes all behave essentially similar to 



P 



AM 



Jm| 



cos (j\m 2 + m(2n + 1)| In f ^~ ll ' 



where the scale A is set by the function \i. Numerical inte- 
gration shows (when /x is majorized by a constant) that the 
numerator scales like the logarithm of the volume while 
the denominator scales as the square root of the volume. 
This is as anticipated due to the explicit factor of l/p 2 - 

Thus, independent of the selected normalization con- 
stants, the second term in F vanishes like ln(U)/^ 1 / 2 , 



leaving again only F(uj 



Since for vortices with 



uj mn Vtu) 2 mn 



fluxes higher than one the level-density does no longer 
vanish at zero, the condition ((31 \ is fulfilled. As these are 
Coulomb gauge configurations, this necessary confinement 
criterion is met. 

It should be noted that the result is not the same as 
on a finite lattice in a thin-center-vortex-only configura- 
tion which seems to indicate a finite value of F(0). 
However, this depends significantly on the volume, and 
from the available ones, it cannot be excluded that F(0) 
does not vanish in the infinite-volume limit also in lattice 
calculations. 



6 Discussion and conclusions 

Summarizing, the spectrum of the Faddeev-Popov opera- 
tor in a one-instanton background field has the following 
features: It is positive semi-definite. It has a continuous 
spectrum of positive eigenvalues, just as in the vacuum, 
but with a different angular structure. The (infinite) de- 
generacy of non-zero eigenvalues does not depend on the 
eigenvalue, and there is no enhancement of eigenmodes 
near zero, but only at zero. Three new zero-modes exist 
at "angular momentum" 1 = 1/2 and 1 = 1. This number 
is independent of the instanton size A, which is the only 
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characteristic quantity of an instanton. Thus the instan- 
ton belongs to the first Gribov horizon, the boundary of 
the first Gribov region. 

Furthermore, the configuration is due to (0 also part 
of the fundamental modular region, and lies on the com- 
mon boundary of the first Gribov horizon and the fun- 
damental modular region. Hence, it is contained within 
the region of field configuration space which is responsible 
for gluon confinement according to the Gribov-Zwanziger 
scenario. Still, the one instanton field does not lead to a 
significant enhancement of the number of zero-modes, as 
required by the Gribov-Zwanziger scenario. It is therefore 
likely not alone, if at all, responsible for or even involved 
in gluon confinement. 

Still, it is tempting to speculate whether an n-instanton 
configuration, consisting of far separated instantons, could 
be able to support n independent zero- modes. The radial 
structure supports this, as at sufficiently far distances ra- 
dial eigenmodes behave like the vacuum solution or vanish. 
But the intricate angular variations, which still exist at 
infinity, makes it far from obvious that a multi-instanton 
configuration does sustain more zero-modes. Only the / = 
1/2 solution shows any sign of localization, and may thus 
be a candidate for a relevant multi-instanton solution. 

Hence, it cannot yet be completely concluded that in- 
stantons are irrelevant to gluon confinement. Therefore, it 
would be interesting to study this connection also in lattice 
gauge theory. This is, however, made complicated due to 
the appearance of sub-structure of topological defects on 
toroidal manifolds, which has been seen e. g. in connection 
with Kraan-van Baal calorons 03] m lattice calculations 
|44|. Digressing, this rises the very interesting question, if 
in these cases any influence on the eigenspectrum of the 
Faddeev-Popov operator will be due to the constituents 
or will arise as a feature of the bound states. 

On the other hand, a vortex of sufficiently high flux 
can sustain an arbitrary number of zero-modes. However, 
a vortex of flux one is not sufficient for this; a larger flux 
is necessary. Only flux three or higher vortices belong to 
the Gribov horizon, while vortices with lower flux are only 
part of the interior of the first Gribov region. Thus, vortex 
configurations, if they belong to the fundamental modu- 
lar region 9 , could contribute to the enhancement of the 
spectrum of the Faddeev-Popov operator, and even be its 
main source. This is again only an enhancement at eigen- 
value zero, as the (infinite) degeneracy of non-zero modes 
is again not dependent on the eigenvalue. 

The fact that a single vortex can sustain a large num- 
ber of zero-modes is different from the instanton case, 
where a single instanton was not able to support more 
than three zero-modes, irrespective of its size. This is to 
some extent surprising, as it has been shown that a vor- 
tex does not carry topological charge when it is oriented 
0, as it is the case in the present calculation. The only 
relevant property seems to be a sufficiently large flux. 

It is noteworthy that the zero-modes, with the excep- 
tion of the 1 = 1/2 one in the instanton case, are not local- 

9 As vortex configurations appear as physical excitations in 
lattice calculations |45j . this seems likely. 



ized, i. e. do not vanish at spatial infinity. Especially the 
vortex zero-modes do still vary appreciably even in radial 
direction at spatial infinity. This indicates that eventu- 
ally non-localized modes may be important in the confine- 
ment problem. This is supported by similar observations 
in recent lattice investigations of the covariant Laplacian 
[4"r>| and also of the Faddeev-Popov operator j^H] . Also in 
the context of the Gribov-Zwanziger scenario, the non- 
localization of such modes play a role . Still this issue 
is not finally settled, and it is not yet really clear, what 
the significance of localization in this context is. 

In addition, the vortex has been found to satisfy one 
necessary criterion for quark confinement in Coulomb gauge 
Thus, if all of these indications are correct, and such vor- 
tices play a significant role in the dynamics of Yang-Mills 
theory and eventually of QCD, they are perhaps the link 
searched for in this work. Nonetheless, it is likely that 
vortices are still only a symptom rather than the origin of 
confinement, as confinement exists also in Yang-Mills the- 
ories with centerless gauge groups like G2 In these 
cases, possibly different types of topological excitations 
play the same role as vortices in SU(2). 

Concluding, the following results have been found: In- 
stantons belong to the common boundary of the first Gri- 
bov region and the fundamental modular region, but only 
support a small number of zero- modes, independent of 
the size of the instanton. This shows that not all config- 
urations on the common boundary provide a significant 
enhancement of the spectrum of the Faddeev-Popov oper- 
ator at eigenvalue zero. Instantons seem again not directly 
relevant to confinement, but it is not clear whether an en- 
semble of instantons can sustain more zero-modes. 

Vortices of flux zero and also of flux one do not belong 
to the common boundary nor even to the Gribov horizon 
and therefore do not play any role in the Gribov-Zwanziger 
scenario. This implies that not all topological field config- 
urations lie on the boundary nor provide zero-modes in 
the spectrum of the Faddeev-Popov operator. Vortices of 
flux 2n + 1 > 3 belong to the Gribov horizon, and lie al- 
ways in the first Gribov region, although pointwise \A^\ is 
not bound for large n. (Oriented) vortices can support An 
zero-modes, and thus provide a significant enhancement 
of the spectrum at zero of the Faddeev-Popov operator. 
This corresponds to results in lattice calculation that the 
enhancement of the spectrum of the Faddeev-Popov op- 
erator vanishes when the vortex content is removed |35j . 
Hints for a related modification of the propagator have 
been found on the lattice .'51 . This indicates that perhaps 
the ghosts are the place to look for traces of topological 
excitations in functional methods. 

Finally, these results demonstrate that certain types 
of topological configuration influence the spectrum of the 
Faddeev-Popov operator, and can give rise to the appear- 
ance of zero-modes. Furthermore, it has been shown that 
this strongly depends on the type of configuration, and 
that even some configurations do not provide an enhance- 
ment at eigenvalue zero. These results provide a first step 
towards an analytical understanding of the intimate link 
apparently existing between a confinement dominated by 
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topological aspects and the Gribov-Zwanziger scenario, 
which are both by now well supported by various calcula- 
tions. In addition, the results presented here also support 
similar findings in numerical lattice calculations. 
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